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Abstract: The D3/D7 system hole-graphically describes an J\f=2 gauge theory which 
spontaneously breaks a chiral symmetry by the formation of a quark condensate in the 
presence of a magnetic field. At finite temperature it displays a first order phase transition. 
We study out of equilibrium dynamics associated with this transition by placing probe D7 
branes in a geometry describing a boost-invariant expanding or contracting plasma. We 
use an adiabatic approximation to track the evolution of the quark condensate in a heated 
system and reproduce the phase structure expected from equilibrium dynamics. We then 
study solutions of the full partial differential equation that describes the evolution of out 
of equilibrium configurations to provide a complete description of the phase transition 
including describing aspects of bubble formation. 
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1. Introduction 

Thermal phase transitions are a crucial aspect of the evolution of the Universe after the Big 
Bang and also in the physics of heavy ion collisions. We have traditionally lacked tools to 
study these transitions though in strongly coupled systems such as QCD. The AdS/CFT 
correspondence |1J, which gives a weakly coupled string/gravity description of a class of 
strongly coupled gauge theories, offers the chance to study similar transitions in detail. 

In this paper we study a simple dual of a theory with gauge fields and quarks which 
has (in the presence of a magnetic field) a first order chiral symmetry restoring phase 
transition as it is heated. Previous analysis of this gauge theory (using probe branes 
in an AdS-Schwarzschild black hole geometry) has explored the first order transition for 
equilibrium (time independent) configurations. The heating or cooling of the system can 
be studied thanks to the boost-invariant expanding or contracting plasma geometry of 
Janik and Peschanski |^]. That geometry, which we will review below, has a moving black 
hole horizon describing the changing temperature in the gauge theory. We will study 
probe branes in this geometry to learn more about the first order phase transition out of 
equilibrium. 

The particular duality we concentrate on is the simplest example of holography with 
fundamental quark fields ||, ^ ^, Q, gj. We do not consider the specific degrees of 
freedom of the theory too crucial - it is some strongly coupled gauge theory that displays 
a chiral phase transition. We hope, in the spirit of AdS/QCD models |9|, |l0fl , that it 
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reflects broad aspects of many strongly coupled systems. The specific gauge theory is 
constructed from the D3/D7 system in type IIB string theory which we will describe further 
below. The theory is the large N c J\f = 4 U(N C ) gauge theory with a small number of 
quark hypermultiplets. We will work in the quenched approximation ]|] (appropriate when 
Nf <C N c ) which on the gravity dual side corresponds to treating the D7 branes as probes 
in the metric generated by the D3 branes. There is a U(l) chiral symmetry (a remnant of 
the SU(4:) R-symmetry of the M = 4 theory) which is broken when a quark condensate 
forms 0, [ll|, |lj]. Several mechanisms for triggering this condensation have been explored. 
The cleanest is when a background magnetic field is introduced [13, 14, 15, 16, 17 1 so we 
will use that mechanism here. The reader might want to loosely view the B field as simply 
the introduction of a conformal symmetry and supersymmetry breaking parameter that 
triggers the strong dynamics to cause the symmetry breaking. Physically, magnetic fields 
may be strong during structure formation in the early universe, in particular during the 
epoch of the QCD primordial phase transition, and in non-central high-energy heavy ion 
collisions studied at RHIC. Running of the coupling in the holographic dual also causes 
quark condensation as has been shown in back-reacted dilaton flow geometries |l8| and 
models with a phenomenologically imposed dilaton profile [19]. The quark condensate can 
be determined in these models and an effective IR quark mass is generated. The theories 
display a massless pion-like Goldstone field and a massive sigma field (since we are at large 
N c it is stable) that is the effective Higgs particle. 

The equilibrium finite temperature behaviour of the theory with a magnetic field has 
been studied in |jl3"l , |i"7)[ . Finite temperature can be included through the presence of an 
AdS Schwarzschild black hole. At a critical temperature the D7 embedding flips from a 
chiral symmetry breaking embedding away from the horizon to a symmetry preserving 
embedding that enters the horizon. The transition is therefore also associated with meson 
melting p(]] - for embeddings away from the horizon there are regular linearized fluctuations 
describing the meson spectrum of the theory. For the embedding that enters the horizon 
there are only in-falling quasi-normal modes describing unstable plasma fluctuations. In 
terms of the quarks of the theory the high temperature phase is analogous to the quark 
gluon plasma phase in QCD whilst the low energy phase is more akin to the hadronic phase 
of QCD. It should be noted though that the gluonic degrees of freedom deconfine at any 
finite temperature so the analogy is imperfect. 

The crucial extra ingredient we shall add to this story in this paper is provided by 
the boost-invariant expanding or contracting plasma geometry of Janik and Peschanski |2j . 
This geometry has a black hole whose horizon moves away from the boundary in time as 
the N = 4 plasma it describes expands and cools. The time reversed solution at zero 
viscosity describes a heating plasma and we will find it useful to discuss that scenario too 
below. The geometry is a late time expansion (when the black hole is small) in powers of 
inverse time. However, by controlling the strength of the magnetic field on the D7 probes, 
felt by the quarks, we can arrange to place the chiral phase transition at any point in the 
evolution so the expansion is sufficient to fully study the transition. 

We will place a D7 brane 1 into the expanding plasma geometry and determine the 
By placing a fundamental string, which corresponds to a heavy quark, into the expanding plasma 
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partial differential equation (PDE) that describes the time dependence of its embedding. 
A good first approximation to the transition behaviour is provided by the equilibrium 
results with the temperature replaced by the temperature as a function of time from the 
moving background. In fact we will see that that is an extremely good approximation 
when talking about the slow or adiabatic heating or cooling of vacuum configurations at 
temperatures even close to the phase transition. If the chiral symmetry breaking embedding 
is heated though, the local minima in the effective potential associated with that embedding 
is eventually lost and the configuration becomes an out of equilibrium configuration. The 
equilibrium results can not describe the subsequent evolution. Similarly excited vacuum 
states must be studied through the full PDE system. 

We first turn to an approximation to the PDE. The solution can be power expanded in 
inverse powers of time. This reduces the PDE to a system of ordinary differential equations 
(ODEs) that are much easier to solve. This analysis was first done in ^jj where the time 
evolution of the high temperature phase of the pure M = 2 theory was studied. Using this 
technique we solve the ODE system and find the moving D7 solution. This method assumes 
that the initial vacuum state is exactly a maximum or minimum of the zero temperature 
effective potential. Again the heating is essentially adiabatic in nature. The expansion 
breaks down if the super-heated state ceases to be an extrema of the effective potential. 
This method allows us to confirm the success of the equilibrium derived results although 
in fact the expansion breaks down before the equilibrium results deviate from the full PDE 
solutions. 

The most interesting out of equilibrium questions lie beyond the adiabatic approxima- 
tion though. In a physical first order transition quantities such as the condensate do not 
jump but the vacuum state instead performs a fast roll from one configuration to another. 
The timing of that transition can and most likely will be spatially dependent i.e. bubbles 
of the true vacuum will form and grow. We turn to solving the full PDEs to study these 
phenomena (with care this can be done using in built PDE solvers in for example Math- 
ematica). For example we are able to watch the super- heated chiral symmetry breaking 
phase roll to the symmetric phase. We can also simulate initial conditions with some extra 
energy (which might for example come from thermal fluctuations) and see such configu- 
rations transition away from the super-heated vacuum before the meta-stable vacuum has 
disappeared in the effective potential. This allows us to confirm some elements of the tran- 
sition such as the length of time in which there is a mixed phase. We also are effectively 
watching a very large bubble form. 

Our main result then is to have developed numerical techniques that let us reproduce 
the phase structure using the PDE solutions and to describe non-equilibrium physics that 
is necessarily present in the first order transition. 

In our final section we also analyze the ODE expansion approximation to the PDE 
solutions for the D7 embedding to make clear the full dependence of the solutions on the 
magnetic field value. In that case the dependence is available analytically. We also show 
the effects of the viscosity. 



geometry, the diffusion constant |2l| and drag force [E2j was computed. 
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2. Holographic Descriptions 



In this section we review the J\f = 2 gauge theory with a magnetic field and its holographic 
description. We discuss the theory's finite temperature chiral phase transition in an equilib- 
rium description. We then review how to study flavour physics in a nonequilibrium set-up 
using holography. 

2.1 Chiral symmetry transition 

We will study the D3/D7 brane model at finite temperature and with a magnetic field fl3| , 



17]. The magnetic field, which causes chiral symmetry breaking, competes with the tem- 
perature that prefers to restore chiral symmetry. Consequently there is a first order phase 
transition. 

The N = 4 SU(N) gauge theory at finite temperature has a holographic description 
in terms of an AdSs black hole geometry which can be written as 

2 r>2 

ds 2 = -^{-g t dt 2 + g x dx 2 ) + -?(dp 2 + p 2 dn\ + dL 2 + L 2 d^ 2 ) , (2.1) 
R z w z 

where eft is a U(l) angle and w = \/ p 2 + L 2 , p := ws'mcp, L := w cos cj) and 

/ 4 4 "i2 4 i 4 

_ (u> 4 - w%y _ w +w H 

Note R 4 = 4:7rg s Na' 2 , and the temperature is given by wh '■= irR 2 T. 

Quenched (Nf <C N) M=2 quark superfields can be included in the M = 4 SU(N) 
gauge theory through probe D7 branes in the geometry The D7 probe can be described 
by its Dirac-Born-Infeld (DBI) action 

Sdbi = ~T D7 J d 8 ^-det(P[G] ab + 2Tia'F ab ) , (2.3) 

where P[G] a b is the pullback of the metric and F a b is the gauge field living on the D7 world 



volume. We will use F a b to introduce a constant magnetic field [ 13 1 



F12 = -F 2 i = B/(W) . (2.4) 

We embed the D7 brane assuming only p dependence: L{p) at constant (j). The full DBI 
action we will consider is then 

Sdbi = J dfC{p) = (J e 3 J dtdx^j J dp C(p) , (2.5) 

where e 3 is a volume element on the 3-sphere and 



_ £ = n = Nf T D7 (RVB) 4 f (l - ^(l + L'^ 



l+^X + i,} (2.G) 
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Figure 1: The equilibrium description of the first order chiral phase transition in the Af = 2 
gauge theory with a magnetic held: (a) The D7 brane embedding profiles L(p) as a function of 
temperature wh- The low temperature symmetry breaking embedding is shown on the left whilst 
the symmetric flat embedding on the right is preferred at high temperature, (b) A plot of the chiral 
condensate as a function of the temperature with the dotted line indicating the position of the first 
order transition. 



with the dimensionless variables defined as 

@Xd:=(^,-±=,-e=) • (2-7) 
\rV2B RV2B RV2BJ 

Note that we can use the magnetic field value as the intrinsic scale of conformal symmetry 
breaking in the theory - that is we can rescale L and p by B. The Euclidean on-shell 
Lagrangian (— C) is interpreted as the free energy density (f2). 

In all cases the embeddings become flat at large p taking the form 

~ c „ 2ira'm q „ (2W) 3 

Lip) ~ m + , m = ==.- , c = (qq) == — , (2.8) 

y " P 2 RV2B W '{RV2B) 3 

where fh and c are identified with the quark mass and the quark condensate, respectively. 
Since we are interested in spontaneous symmetry breaking we impose fh = in the UV 
(p — > oo). Then the value of c is determined by requiring regularity in the IR {p 1). If 
there is more than one solution then we choose the one giving the minimum free energy. 

The results are displayed in Fig. [l]. At low temperatures T <C B, the black hole is 
small and the embeddings are repelled from the origin of the L — p plane (Fig. [I] (left)). 
This behaviour is a result of the inverse powers of w 4 , when wh <C 1, in the last term 
of the Lagrangian (|2,6|), which causes the action to grow if the D7 approaches the origin. 
Consequently c is non-zero and chiral symmetry is broken. If the temperature is allowed 
to rise sufficiently then the black hole horizon grows to mask the area of the plane in which 
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the inverse w 4 terms in the Lagrangian are large. At a critical value of T (% = 0.2516, 
Fig. H (middle)) the benefit to the m = embedding of curving off the axis becomes 
disfavoured and it instead lies along the p axis - chiral symmetry breaking switches off. This 
transition is first order since the condensate vanishes discontinuously, which corresponds 
to the embedding change from the blue to the red one. At higher temperatures (Fig. [I] 
(right)), the embedding stays flat, c = 0, as expected in the chirally symmetric phase. The 
corresponding condensate vs. temperature (— c-wh) plot is shown in Fig. [[] (bottom). 

2.2 A boost-invariant expanding plasma 

The geometry (|2.1| ) is dual to a system in thermodynamical equilibrium and therefore not 
suitable for the description of the chiral phase transition in a rapidly expanding plasma, 
in which the transition is basically a non-equilibrium process. Boost-invariant expanding 
M = 4 SU (N) plasmas out of equilibrium can however be described by the time-dependent 
background found in Q (see [^4|] for a review) . In the following we review the basic features 
of this background and discuss the embedding of probe D7 branes dual to quenched flavours 
("quarks") in the plasma. 



The boost-invariant geometry is a 5D spacetime with coordinates {r, y, x±(= 



x , x 



which, apart from the holographic direction z, parameterise the 4D spacetime on the bound- 
ary. The longitudinal position plane is parameterised by the proper time r and rapidity 
y (related to x ' 3 as x° = rcoshy, x 3 = rsinhy), the transverse coordinates are collected 
in x±. We also add a five-sphere to obtain a full type IIB supergravity background. The 
metric is then of the form 



1 



_ e a(r,z) dT 2 + e b(r,z) T 2 dy 2 + e <r,z) ^2 



dz 2 



2 

5 > 



(2.9) 



R 2 z* \ -/ 
where R is the radius of the AdSs space. At late times, the coefficients can be expanded 
to first order as 0, 63, E26| 
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(2.10) 



r 4/3 



1 



U - T l/3°0 > 2 l/2 3 3/4 

eo is a free parameter of mass dimension 8/3 and is related to the energy density, while 
770 is related to the shear viscosity, v is a scaling parameter valid at large r. Note that 

00 in powers of 1/r 2 / 3 as § |2E 



(2.11) 



a(r, z),b(r, z) and c(t,z) are expanded around r 

00 / 

a(r,z) = 



1 



a n . (v 



n=0 



.3/8 



26 1 
(2.12) 
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and similarly for b(r, z) and c(t,z). The coefficients are functions of a scaling parameter 
v only. Due to this scaling behaviour the complicated PDE Einstein equations can be 
reduced to an ODE system, which allows analytic solutions. At early time r <C 1 there is 
no scaling behaviour and we should solve the full PDEs. 

To find which gauge theory state corresponds to this bulk metric according to the 
AdS/CFT dictionary let us expand the metric around z = 0, 



9yy 



a a(r,z) 



T 2 e b(r, Z ) 



2vr 2 . 

c 



2rj 

T 



+ 0(z 6 ) 
2rj 



(2.13) 



where 2 



9ii(= 922) ^e c ^ 



N 2 C e 



1+ Nf 



+ 0(z 6 ) 



2tt 2 r 4 /3 



2tt 2 



3/4 



(2.15) 



The leading terms (of order O(z )) of the metric elements ( 2,13| ) is simply the Minkowski 
metric in the r— y coordinate system, and the subleading terms (of order 0(z )) correspond 
to the expectation value of the energy-momentum tensor, i.e. 



00/ 



^ IT ) 



22/ 



3 



(2.16) 



This energy-momentum tensor is precisely that of a longitudinal viscous boost-invariant 
N = 4 SYM conformal plasma with finite rj. 

In g TT we recover a time-dependent emblackening factor, which describes a moving 
horizon. The size of the horizon determines the time-dependence of the 'temperature' 
as p| pi 



T(r) 



V2_ 

R 2 



TT 



f) 



1/4 



7TT 



1/3 



24 /4 r 2 /3^ 



(2.17) 



If we assume that the time-dependent entropy density s(t) has the same form as in the 
static case, a(r) = (vr 2 /2)A^ c 2 T(r) 3 , then the ratio rj/s can be computed as |26| 



1 
s 



4tt 



+ 0(t 



2/3^ 



(2.18) 



which at large r agrees with the known static bound. Note that the numerical value of rjo 
in ( 2.11 ) is crucial to get l/4ir. 



j^y- has been factored out to apply the AdS/CFT dictionary 



N 1 



(2.14) 
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Finally we note the uniqueness of the gravity solution. There is a potential singularity 
at v = 3 1 / 4 in (|2TT0|) . The curvature invariant, R^^R^P' 7 = R (v) + R^v)^ 2 ^ + 



R2{v)t 3 / 4 + • • • is only regular at each order if one makes the specific choices above. We 
must choose —1/3 for the power of r in ( 2.11| ) to make Rq regular We must also choose 



( |2.10[ ) with the specific choice of the numerical value 770 as in (2.11) to make R2 regular |^6| 
is always regular and does not give any constraint.) 



Flavours in this background were first studied in J23| . Knowing the explicit background 
geometry ( |2.9| ) we can study flavour physics using the D7 brane DBI action. The action 
schematically reads 



S ~ J d 8 ^-det(P[G} ab ) ~ J drd P Tp 3 AJl + {d p Lf , (2.19) 



where A and IB are complicated but known functions of r and p, see ( |3.5| ) and ( |3.6| ) in the 
next section. L = L(t, p) is the embedding profile of the D7 brane and is assumed to be a 
function of r and p. 

The equation of motion coming from ( [2.19 ) is a non- linear PDE. However it can be 



semi-analytically solved by the late-time expansion 

00 

L(T,p) = m + J2fi(P>~* ■ ( 2 - 2 °) 

i=i 

This expansion inherently assumes that the late time configuration is precisely the equi- 
librium vacuum and when we heat it there is no excess energy. This leads to an adiabatic 
approximation. The fraction 1/3 is chosen because all exponents of r in A and IB are integer 
multiples of 1/3 and the constant m reflects the fact that the embeddings get flat at large 
r. Here we only consider a large bare quark mass (m 3> 1) and Minkowski embeddings, 
i.e. D7 configurations which end well above the black hole horizon. 

This reduces the partial differential equation to an infinite set of ordinary differential 
equations given by 

P -% (p 3 fm = 9(m 8 r + £ ; 2)5 /, , (2.21) 

Is = 1 , J10 = -4t? £o 1/4 , otherwise Ji<i 3 = , (2.22) 

where we do not consider terms with % > 13 since they are beyond the validity regime of 
our approximation of the boost-invariant metric. The asymptotic solution at large p is 

fi ~ mi + ^ • (2-23) 
p 2 

We set mi = so the bare quark mass is zero and not time-dependent. We also impose the 
condition //(O) = for regularity. With these boundary conditions, we find fi = except 
for /g, fiQ. As a result, one gets 

L(T,p) = m + c r r , (2.24) 

(m z + p^Y 
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where 

c(r) = ( r -8/3 _ 4r?oe -V V lo/3 + •••). (2.25) 

The condensate approaches zero as ~ — r™ 8 / 3 and the viscosity has a "dragging" effect 

^ + r -10/3. 

We close this section with some technical remarks on black hole embeddings. In the 
Feffer man- Graham (FG) coordinate system ( |2.9| ) one cannot approach the horizon for fixed 
r as in the Schwarzschild black hole metric. There is also an extended background written 
in terms of Eddington-Finkelstein (EF) coordinates, where the spacetime is well defined 



across the horizon [27, 28] and black hole embeddings may be described more consistently 
there. However we have found that the embedding configuration is not easy to handle 
in those coordinates. Below, to enable us to study black hole embeddings, we use FG 
coordinates but with a cut-off slightly above the horizon. Since the FG coordinate system 
is a good patch near the horizon at large r, by restricting ourselves to large r, we may 
capture the essential physics of the embedding. This will allow us to go beyond the results 



in 23 . 



3. Out of equilibrium description of the chiral phase transition 

In the previous section we reviewed the holographic description of dynamical flavours in 
an expanding plasma p3j| . A boost-invariant background with embedded D7 branes is 
however not sufficient to study the chiral transition. In that case the D7 embeddings are 
always flat in the chiral limit (m —¥ 0) and the system is always in the chiral symmetric 
phase with vanishing quark condensate. In order to describe the transition to the chiral 
broken phase, we need a repulsive effect to compete against the attractive force of the black 



hole. As in the static case reviewed in section 2T , this can be achieved by turning on a 
magnetic field. In this section, we will therefore consider D7 branes with a world-volume 
magnetic field in the dual geometry of an expanding plasma. In this way we will find the 
holographic dual of a chiral transition and deduce the dynamic effective potential for the 
time-dependent quark condensate. 

3.1 D7 flavour brane action 

The background metric for a boost-invariant expanding plasma can be written as 

ds 2 = ^ { _ e a(r,r) dT 2 + ^(^2^2 + e <r,r) ^ + ^^2 + ^2 + ^2 + L 2 #2) ^ 



where r 2 = p 2 + L 2 . The S 5 part is written as in ( |2.1| ) and the AdSs part follows from fl2.G| ) 
with a change z — > R 2 /r, i.e. a(r, r) = a(r,z — > R 2 /r) in ( |2.10| ) and similarly for 6(r, r) 
and c(t, r). 

We are interested in time-dependent D7 brane embeddings of the type L = L(r, p). 
The corresponding DBI action is 

Sdbi = ~T D7 J d 8 ^-det(P[G] ab + 27ra'F ab ) , (3.2) 



-9- 



where we turn on a constant magnetic field [13|, 



F 



12 



-Fzi = B/(2ira') 



(3.3) 



in order to induce the chiral symmetry breaking. 
More explicitly, the D7-brane action reads 



S = n drdprp 3 A 



1 + 
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3 + tj 4 



-2/3 



Goth" 1 4 T~ 2 / 3 



e K 



T 1 /*^p 2 +L 2 



N = iV/T D7 y e 3 y dyd 2 x ± , 



(3.4) 

(3.5) 
(3.6) 

(3.7) 
(3.8) 



where A 



= a/2+6/2+c 



e~ a , C = e- 2c , r/o = l/(2 1 / 2 3 3 / 4 ) as in ( pUP , and e 3 is the 



volume form of the three-sphere. For vanishing 5, the action reduces to that in ( [2.19D , 
It turns out to be convenient to work with the rescaled variables 



R 



P 



yfBRp , L = VBRL , t = —=t , e = S 4 / 3 R~ 8 / 3 e 



(3.9) 



For R = 1 the action then reads 



S = nB drdprfA 



1 + C- 



1 



L 



1 + L' 2 - I 
(pt + L 2 ) 2 J \ (p^ + L 2 ) 2 

with A, B, C as in (|3.5D-({37^) but now expressed in terms of r, eo,p, and L. 



(3.10) 



3.2 Naive equilibrium based approximation 

An obvious first approximation to understanding the time dependent chiral phase transition 
in this set up is to use the equilibrium results from section 2.1. There we described how 
the D7 embedding behaved in the background of a fixed size black hole and determined 
the quark condensate as a function of temperature. In the cooling plasma geometry of 
section 2.2 the black hole horizon moves as a function of time r as 



a/4 



r H 



31/4^1/3 



1 



2e] /4 f 2 /\ 



(3.11) 



If this heating were very slow (as it is at large r) we would expect to be able to plot the 
quark condensate against r by simply substituting for the temperature T in the equilibrium 
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V(c) 




(a) (b) 



Figure 2: A summary of our main results showing the first order phase transition in our vari- 
ous approaches, (a) The condensate c(r) in the equilibrium (black), adiabatic (dashed) and non- 
equilibrium (blue) approaches. In the equilibrium description, the different branches correspond 
to the extrema of the potential V(c) (b). The complete potential is obtained by a U(l) rotation 
around the vertical axis. The dashed red curve schematically shows the path of the adiabatic 
heating evolution. 

results. We show that plot in Fig. ||a (black solid curve) which follows directly from the 
c-T plot in Fig. [j]. 

We can recast the c vs. r plot schematically as a time dependent effective potential by 
fitting it at each value of r to a potential of the form 

V = m 2 !? + Ac 4 + A'c 6 (3.12) 

using the values of c at the extrema to fix the parameters (the overall scale is not set but 
the figure is intended to only be schematic). We plot this in Fig. Qb. This is the standard 
picture of a first order transition. 

It is important to interpret these plots correctly. Firstly note that there is a c = 
solution for all r. If we begin at high temperature with the symmetric flat embedding and 
cool slowly or adiabatically {i.e. move to large r) then we can remain in that embedding for 
all r. Above the time r a = 20.6 two extra solutions for the condensate develop and these, 
as we will discuss, trigger the first order transition. At large r the c = flat embedding 
becomes a local maximum of the effective potential. To stay in the flat embedding for all r 
is the extreme limit of super-cooling the high temperature phase into an unstable vacuum 
state. 

The top trajectory in the c vs. r plot is best thought about in the time reversed 
solution that is heating up. At early times (large r in the plot) the solution is the T = 
symmetry breaking D7 embedding vacuum. Now as we heat adiabatically the condensate 
tracks along the trajectory to smaller r. The solution ceases to exist at r a , where the cusp 
in the (black) curve is, indicating that the minimum of the effective potential corresponding 
to this solution has ceased to exist at r . At this point the brane will move quickly, as 
an out of equilibrium configuration, ending as the flat embedding or oscillating about 
it. The real-time evolution of the configuration is represented by the path of red dots in 
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Fig. 2b. Starting from the global minimum with c 7^ 0, the minimum is lifted up and the 
configuration moves adiabatically. At r a the minimum disappears and the system roles 
down the potential to the true vacuum with c = 0. The equilibrium results can not tell us 
about this motion. 

Now we also see the role of the middle arc in the plot - these solutions are symmetry 
breaking embeddings that end on the horizon and they are a local maximum of the po- 
tential lying between the two minima (the embeddings of the top trajectory and the flat 
embedding). 

In fact the first order transition point between the symmetry preserving and the sym- 
metry breaking embeddings can be computed from the free energy in the equilibrium 
computation. This transition occurs at r c = 27.5, where the local and the global minimum 
in the effective potential interchange their roles, i.e. for r c > 27.5 the symmetry breaking 
vacuum (with c / 0) becomes the global minimum. Another important time determined 
by the equilibrium computation is % = 34.9 where the mixed phase ends at large r. 

Our first task in solving the PDEs for the D7 brane motion that result from ( |3.4| ) is 
to demonstrate that this description is essentially correct - we will see that it is. We will 
try to find confirmation of the times r a _ c of these events. The more interesting task is 
that we will be able to follow the evolution of a particular initial condition through the 
phase transition. Of course in the first order transition the brane configuration does not 
discontinuously leap between the symmetry breaking embedding to the flat embedding but 
evolves continuously. We will provide solutions for this evolution. 

Another interesting phenomena associated with a first order transition is bubble for- 
mation. In real systems thermal energy will lead to volumes of space which have more 
energy than an equilibrium like state in a local minimum of the potential at any particular 
time. These volumes may "climb" over the potential hill to the other local minimum during 
the mixed phase period shown in Fig. ^. These bubbles then grow or contract triggering 
the phase transition around r c ending any super-heated or cooled phase. We will not look 
at (#3) spatial dependent brane embeddings. However, we can inject kinetic energy in the 
holographic directions of our description into the brane configuration before we heat or 
cool it to see the configurations moving more quickly between the two local minima than 
the lowest energy configuration. This will allow us to test the time period in which the 
mixed phase exists in the out of equilibrium problem. 



3.3 Adiabatic dynamic D7 brane embeddings 

Our first study of the PDEs describing the D7 embedding in the expanding plasma ge- 
ometry will be to study adiabatic expansion. The non-linear partial differential equations 
resulting from ( 3.10| ) can again be transformed into a set of second-order ordinary differ- 



ential equations as we reviewed for the case with no symmetry breaking in section 2.2. For 
that, we use the late-time expansion 

00 

L(r,p) = Mp) + ^2f l (p)r-i. (3.13) 
i=l 
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p 



(a) (b) 

Figure 3: Solutions of the ODEs in the adiabatic approximation - /o,4,s,io(p) f° r ?o — 1- ( a ) /o 
(red) compared to the dashed curve which is the profile obtained in the static case at T — 0. (b) 
/o (red), / 4 (orange), f 8 (yellow), f w (green). 



Thus we will be following the evolution of the symmetry breaking embedding as r decreases. 
There is also the solution fi = which corresponds to the symmetric embedding. Note 
that, contrary to ( |2.20| ), a nontrivial asymptotic embedding fo(p) is assumed because of the 
repelling effect of the magnetic field. The equations for every fi(p) can be obtained order 
by order in r~3, which will be solved recursively. Unlike ( p. 21 ) these equations are quite 



lengthy and will not be presented here. They do not allow for an analytic solution such as 
(2.24) but can be numerically solved without any difficulty. As in the static case 0, there 



are expected to be two types of brane solutions, Minkowski and black hole embeddings, 
depending on whether the brane ends on or above the horizon. 
For Minkowski embeddings, the boundary conditions are 

#(0) = and /i(oo) = 0, (3.14) 

as before. With these boundary conditions, we find non-trivial profiles /, = fi(p) for 
i = 0,4, 8, 10, and fi = otherwise. A general dynamical embedding function is therefore 
of the form 

L(r,p) = f (p)+ y4/3 + ?8/3 + , (3.15) 

where we neglect terms with i > 11 since they are beyond the validity regime of our 
approximation of the boost-invariant metric. 3 The numerical plots of the non-trivial profiles 
fi (i = 0, 4, 8, 10) are shown in Fig. |[ 

These profiles have the following qualitative properties. As compared to the solution 



( 2.24 ) for B = 0, there are additional non-trivial profiles, /o and f^. The profile /o 
agrees with that in the static case (r — > oo) jO^] (dashed curve in Fig. |3|a) We found 
that the equation for /o is independent of £q and 770, which is natural since at r — > 00 the 
system achieves equilibrium at low temperature and will not depend on the non-equilibrium 
dynamics (r/o) or specific initial conditions (eo) anymore. The profiles f± and f$ depend on 
£q and have a negative sign, which reflects the attraction of the D7 brane to the black hole. 



3 In principle there may be a finite /12, but we will ignore it since it is a higher order term. 
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(c) t = 22 (d) r = 43 



Figure 4: Embedding profiles L(r,p) obtained from the adiabatic ODE expansion method for 
6q = 1, rjo = 0.3102 before (a), during (b,c) and after the phase transition (d): Minkowski embedding 
(green), black hole embedding (red), flat embeddings (blue). 



/io is the first profile which depends on r]Q. In contrast to f± and f$, fio has a positive sign 
showing the dragging effect, as in the zero B case around Q2.25 ). The big amplitudes of 
fi, fs, and /io will be suppressed at large r by negative powers of r. For self-consistency 
we will only consider the r region where all the sub-leading terms are well dominated by 
the leading terms, i.e. fo 3> /4T -3 ' 4 3> • • • . Some of the final embedding profiles L(r, p) 
are shown in Fig. |], where the green lines are plotted by plugging the numerical data of 
Fig. I into ( gig ). 

We would not expect this expansion approach to work beyond the point where the 
symmetry breaking embedding ceases to be even a local minimum because higher order 
terms will grow. One can nevertheless plot solutions using fo — f\o that reach down all 
the way to the horizon. We will use the full PDE solutions in the next section to test the 
point where the expansion has broken down. 

In fact early in our studies we tried to use the expansion even for black hole embeddings. 
The ODEs for the do not contain the horizon however. We attempted to put the horizon 
in by hand by imposing boundary conditions relevant to a black hole. For each point in 
the L — p plane if we assume a horizon is present we can deduce the value of r from ( 3.11 ) 
- one can impose some boundary condition such as orthogonality to the posited horizon 
and seek solutions for each /j. This correctly gives one massless solution at each r and 
they look very similar to the equilibrium black hole embeddings. In fact though on solving 
the full PDEs we realize that this is far beyond the point where the expansion method has 
collapsed! We will include some of the resulting curves below though as evidence of the 
break down of the ODE approximation. 

Fig. H| shows the embedding profiles at various stages of the evolution of the expanding 
plasma. The black hole horizon ( [3.11 ) is indicated by a black quarter circle in the figure. Its 
size decreases with time corresponding to a cooling and expanding plasma. Fig. |I| reflects 
the division of the quark-gluon plasma into the three phases described above. 
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The quark condensate as a function of time can be read form the asymptotic form of 
the solutions 

L(r,j5)~^. (3.16) 
P 

We plot these results in Fig. ^a (the dashed lines) for comparison to the equilibrium inspired 
results. The main result here is that the late time behaviour is indeed just the equilibrium 
expectations. Where the curved line deviates from the equilibrium results is in fact a 
sign that the expansion used in this section has broken down, the configuration is out of 
equilibrium and adiabatic behaviour is no longer possible. Full solutions of the PDEs will 
show this in the next section. 

3.4 Full PDE solutions 

In the previous section we interpreted our solution obtained from the ordinary differential 
system as the evolution of extremum states in the potential. This evolution also reflects 
the real time dynamics of the plasma whenever there is a unique deep global minimum for 
the vacuum state, because the embedding time dynamics is expected to be well localized 
around the minimum. This is the case in the chiral symmetric phase at early times and 
in the broken phase at late times. In the mixed phase at intermediate times there appears 
also a local minimum in addition to the global minimum. The small potential barrier 
separating these minima is potentially easily overcome by a fluctuation. The path of the 
local maximum seen in the adiabatic approach will not be realized as a dynamical solution 
since it is unstable. Further below the critical time r a where the symmetry breaking 
minimum disappears a heating vacuum will be left in a very non-equilibrium state that 
again can not be followed using the expansion technique of section [3.3| . 

To find the real time evolution of the chiral transition out of equilibrium, we need to 
solve the PDE directly and compare the solution with our previous approximate solutions. 
In practice it is more convenient to consider a heating process than a cooling one 4 because 
we can use a well defined starting configuration at large r as an initial condition of our 
partial differential equation, i.e. 

L(t ->■ oc,p) = f , d?L(r -> oo,p) = , (3.17) 



with /o as in ( 3.15j ). For simplicity, we impose only Neumann boundary conditions at p = 



associated to Minkowski embeddings and a zero bare quark mass condition at p — > oo to 
study the spontaneous symmetry breaking, 

dpL(r,p = 0) = , L(r,p ^oo) = 0. (3.18) 



The conditions ( 3.17 ) and Q3.18 ) completely determine the dynamics with the partial dif- 



ferential equation derived by varying the action flOpp . A gain, the actual expression of the 



equation of motion is lengthy and will not be shown here. 



4 In principle, time reversed heating is justified only at zero viscosity since a finite viscosity results in 
decreasing entropy. However, we keep a finite viscosity in our numerics since it has a negligible effect and 
the results are relevant to the cooling case. 
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Figure 5: The embedding profiles L(t, p) from the full PDE solution for a solution starting in 
the low temperature symmetry breaking vacuum (k = 1). 

L[p,r\, rjO = 0.310202 




Figure 6: 3D plot of the embedding profile L(t,p). 

Numerically we solve the equation using Mathematica's inbuilt PDE solvers. These 
are somewhat temperamental and one needs to spend considerable time adjusting precision 
tolerances in order to find smooth solutions in sensible periods of computer time. When 
we have such solutions we test their stability to changes in precision settings to ensure they 
are reliable. 

With the boundary conditions above we can run our simulations until the embeddings 
touch the black hole. Beyond that one needs dynamic boundary conditions along the black 
hole surface. At least in the coordinates we use here this is a hard problem. We have found 
a simple trick that seems to produce sensible black hole embeddings though. After the 
D7 has touched the horizon at p = we artificially hold the embedding at the top of the 
horizon. The large p evolution of the D7 is local and relatively unaffected by this incorrect 
embedding at small p. Further, as has been observed before, solutions shooting from the 
black hole horizon experience a numerical attraction onto the unique regular black hole 
embedding ending at a given point on the horizon. The result is that we get numerical 
solutions like those shown in Fig. [5] (e.g. the red, blue, or orange curve) where the D7 
follows the horizon before shooting out to large p. We believe that these represent very 
good approximations to the large p embeddings solutions. Since we extract the condensate 
c at large p we will live with the improper near horizon behaviour. It would of course be 
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Figure 7: Plots of the IR position of the D7 brane against time for a number of large time initial 
conditions with different energy. Initial velocities: black (k — 1), red (k — 10), green (k = 30), blue 
(k = 40), purple (k = 50). The thick black line is the horizon. 



interesting to try to improve on this with dynamic boundary conditions in the future. 

As a first example of a solution we will study the super-heated symmetry breaking 
vacuum. At large r we use the leading terms in the expansion from Section 3.2 to find the 
UV configuration - one needs to use several terms in the expansion to find the embedding 
with no extra energy. We then solve for the evolution to low r. In Fig. || we show plots of 
the embedding L(p) for different r with the black hole's position for each r also shown. We 
expect that the near horizon behaviour is not correct but the far UV embedding should 
approximate the solution we seek well. In the 3d plot of L(t, p) shown in Fig. ^ this 
corresponds to excising the interior of the region indicated by the dashed blue line. 

A smooth evolution is apparent. To compare this to our various approximations above 
we also plot the condensate as a function of time in Fig. ^a (blue curve). Again the solution 
follows the equilibrium estimate and the expansion solution at large r. In the period T a _ c 
it follows the equilibrium result not the ODE expansion results showing that expansion 
had broken down before the brane left the mixed phase. The success of the equilibrium 
approximation suggests we should take its estimate of the transition points r a (where the 
local symmetry breaking minimum vanishes) and r c (where the two minima of the mixed 
phase are degenerate in energy) as correct. The full PDE solutions allow us to know in 
addition the behaviour of the condensate when we have heated above the temperature 
where the super-heated phase has stopped having a local minima (beyond r a ). This is the 
main result of our analysis here. 

It would be nice to test the equilibrium configurations estimate for the length of time 
in which the mixed phase exists (?& — r a which is about 10 units of t). This we can do 
by looking at some simple out of equilibrium configurations. For the plots so far shown 
we computed L(p) at some large r from the term in the expansion for L(p, r). We can 
give the configuration more energy by simply multiplying that L(p) by a numerical factor, 
k 5 . One might think of these initial states as thermally excited versions of the asymptotic 



J In practice, we fix some r max and modify the initial conditions (3.17) as L(r — > T max , p) = L 



' ODE/ 



and L(j -> r max ,p) = n * I/ ODB (r max , p), where L ODE is given by (|3.15|) 
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Figure 8: The IR speed of the D7 brane against time for the k = 1 initial condition. 



vacuum. In Fig. we plot the evolution of L(p = 0) as a function of t for a number 
of such states. At large r when the theory is cold the unique vacuum is the symmetry 
breaking one and with extra energy the configuration oscillates about that minimum. The 
motion is simple harmonic as can be seen from the independence of the period on the 
amplitude of the oscillations. As the solutions approach r c the different solutions begin 
to diverge. The solution with k = 10 lies close to the equilibrium curve - the oscillations 
about the minimum are small and the state stays super-heated. When k = 30, 40 L(0) 
falls more quickly suggesting that at least some of the brane's length in p has escaped 
the local potential minimum. The upwards wiggles suggest that some of the length is still 
repelled back into the well though. Finally though by k = 50 the brane has ridden over 
the potential barrier and escaped the local minimum. The difference in arrival times at 
the horizon for these configurations (about 6 units of r) is a rough estimate of the period 
of the mixed phase. It seems to broadly match the equilibrium inspired picture again. 

In Fig. H we plot the motion of the IR p = end point of the brane with time for 
a low energy configuration (k ~ 1). The brane certainly seems localized in a minimum 
down to a time of order r = 27 (compare to the equilibrium estimate that the mixed phase 
begins at % = 34.9 and becomes meta-stable at r c = 27.5). Further the steepest period of 
acceleration is below r = 22 (to be compared with the equilibrium estimate for when the 
metastable vacuum ceases to exist r a = 20.6). We conclude that this system is indeed a 
super-heated state that survives in the local minimum until very close to the equilibrium 
estimate for r a . Note the other obvious feature in the plot is the deceleration just below 
r = 20. This corresponds to where in our simulation the D7 first impacts on the black 
hole - at smaller r we hold the D7 at the horizon as discussed above so this behaviour is 
an artefact. A full solution would continue to accelerate along the horizon. 

These configurations with excess L are also very much linked to bubble formation. 
A bubble forms in the mixed phase when a volume of space has excess energy due to a 
thermal fluctuation and escapes the local minima early. Here by treating the whole space 
as one we are essentially describing the formation of a large homogeneous bubble. It would 
be interesting in the future to try to study x dependent initial conditions to understand 




how quickly or slowly bubbles grow. 

4. Dependence of the condensate on B 

In the previous section we focused on the chiral transition induced by a magnetic field, 
which simply played the role of an intrinsic symmetry breaking scale. In this section we 
turn to the physics depending on B more quantitatively. 

To make analytic progress we will concentrate on the adiabatic (ODE) approach in 



which we rescaled all variables by some power of B, see (3.9). In order to study the 
dependence of the quark condensate on B, we use the original parameters, in terms of 
which the condensate can be expanded as 

r(T n\ _ R 3/2 ( r , 1 c 4 (B-^e ) 1 cs(B-^e ) 
c(t,B)-B< c + -^ + 



1 c w {B-We ,Tio) 

\g5/3 7-10/3 



(4.1) 



Note that Co is independent of £o,Vo and r/o enters only in cio- The B dependence of the 
leading term agrees with that in the static (zero temperature) case [13], i.e. it scales as I? 3 / 2 . 



The first subleading term (~ r~ 4 / 3 ) may be compared to the finite temperature case [|l3|, |17|]. 
In the adiabatic approximation, where T ~ t -1 / 3 , this term scales approximately like 
T 4 i?~ 5 / 6 . However, due to the B dependence of C4, this scaling is not exact. In general, 
the effect of the subleading terms is to lower the exponent in c(r, B) ~ B v to a value 
v < 3/2. The scaling of the total condensate c(t,B) with B will again be determined 
numerically. 

Fig. |^a shows our results for the condensate for various values of the magnetic field. 
We have only plotted the condensates associated with the Minkowski embeddings, which 
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Figure 10: r* as a function of B (eq = 1) from the ODE adiabatic approximation results. 



exist for some time r > t*(B). Here t* is denned as the time, when the Minkowski solution 
meets the black hole solution (This corresponds to the cusp in the dashed curve in Fig. §a) . 
We take r* to mark the time of the phase transition. As can be seen from Fig. ||a the early 
time behaviour is an imperfect approximation to the full PDE solutions but the later time 
solutions are consistent with the solutions. We find that for fixed r, or equivalently for fixed 
temperature, the condensate grows with increasing magnetic field. In the limit r — > oo (i.e. 
at zero temperature), we find the following dependence on B: 



lim c(r) = c B 3/2 = 0.2235 3/2 , 



(4.2) 



which is in agreement with the zero temperature result |13[. For earlier times, the depen- 
dence on B is shown in Fig. ^b. We find numerically that c(t, B) ~ B u where the power v 
decreases from 1.5 at large r to approximately 1.0 at small t. In other words, the depen- 
dence on B tends to become linear at high temperatures. Our results hold for sufficiently 
strong magnetic fields. At small B the system is in the symmetric phase (c = 0). 

The tendency for the condensate to increase with B is in qualitative agreement with 
observations in chiral perturbation theory |3(| (oc B 3 ! 2 for strong fields, oc B for weak 
fields), in the Nambu-Jona-Lasinio model |3lJ (oc B 2 ), in a confining deformation of the 
holographic Karch-Katz model gjj (« B 2 ), and in SU(2) || (oc B) and SU(3) @ lattice 
calculations (oc £> 1,6=t0 ' 2 ). The dependence on B typically ranges from linear to quadratic 
behaviour, i.e. the powers of B are in the range 1 < v < 2. 

We may also study the effect of B on the time of the chiral transition which is marked 
by r*. From the D-brane picture we expect that r* decreases with B. For large B, the 
repelling force caused by the B-field is much stronger than the attractive force of the 
black hole. Even at early times, Minkowski embeddings, associated to stable mesons, 
are therefore favoured over black hole embeddings, which implies that the meson melting 
process sets in at some earlier time, i.e. at higher temperatures. 



Fig. 10 shows function of B. We find that r* indeed decreases with increasing B 

as t^B) ~ £?~ 155 . The numerically found exponent —1.55 for the scaling of B is close to 
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Figure 11: The condensate c(r). The green curve is the leading term, the blue is up to subleading 
term, the red is up to third term and the black is the whole expression. 



— 1.5 for the case of vanishing shear viscosity, which can be explained as follows. At r = t* 
and T] = 0, the horizon is located at 

-1/4 1/4 1 
f =^0 = ^0 r "l/3 = I f43) 

r * 31/45:1/3 31/4^ T * 3V4(i2.45)i/3 ' ^ 

where the last equality is from the numerical value of f* at £q = 1. Thus 

r* = 12.455- 1 - 5 , (4.4) 

where for the numerical analysis we chose = 1. The deviation from —1.5 is due to 
the effect of the shear viscosity. We also numerically confirmed that the exponent is 1.5 
without viscosity. 

This scaling of r* may be compared to results for the critical temperature T c of the 
phase transition in the static case. In the adiabatic approximation (when r/o = 0), T ~ 
r~i/ 3 , which implies T* ~ B 1 ^ 2 . This square root behaviour is in agreement with the result 
for the critical temperature T c in the static approach fll7| and with [ 30 1 for strong magnetic 



fields. It is also in qualitative agreement with studies in QCD [34, pq] . 

Finally we consider the effect of changing the viscosity. The viscosity effect is very 
small, since it is doubly suppressed by both large r and small tjq. This can also be seen in 



Fig. 11, where the condensate is plotted for four cases: the green is the leading term, the 
blue includes the subleading term, the red includes up to the third term and the black is 
the full expression. The viscosity effect is then the difference between the red and black 
curves. There is a small visible difference only at small r. 



5. Discussion 



We have analyzed the first order chiral and meson melting phase transition in a warming 
or cooling strongly coupled gauge theory with quarks using the AdS/CFT Correspondence. 
We have developed numerical techniques to study the PDE that describes the motion of 
a D7 brane in a time dependent geometry. In particular this allows us to explicitly find 
smooth solutions of the non-equilibrium configurations that are necessarily part of the 
transition. These results confirm the equilibrium analysis of the transition but also go 
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beyond them. For example we have described the formation of large homogeneous bubbles 
in the mixed phase of the transition. In the future it would be interesting to study spatially 
inhomogeneous bubbles to understand their growth although this would require the solution 
of a 2+1 dimensional PDE which is potentially more numerically tricky. 

To keep our problem a simpler 1+1 dimensional PDE we also restricted motion of the 
D7 brane to the holographic coordinates L and p. There could also of course be fluctuations 
in the holographic angular direction which we have not described. Such configurations 
might be useful in the study, for example, of disordered chiral condensates [j36|. We may 
also consider the finite density case by turning on the time component of the U{\) gauge 
field on the probe brane [17, 37]. 

We hope that the system we have studied can shed some light on first order transitions 
in a range of strongly coupled gauge theories including perhaps QCD. The M = 4 back- 
ground, however, does deconfine in the presence of an infinitesimal temperature, which is 
not expected in simple QCD. We note though that there are several arguments in favour of 
a possible splitting of the deconfinement and chiral transitions in QCD in the presence of a 
strong magnetic field 34]. While the temperature of the chiral transition increases [38, |3~4], 
the temperature of the deconfining transition either decreases with increasing magnetic 
field [39] or increases but much slower than the chiral transition [34]. Both transitions 
become first-order transitions and a new phase with broken chiral symmetry and decon- 
finement appears for sufficiently strong magnetic fields. So, it is the region of the T-B 
phase diagram with strong magnetic fields and above the deconfinement temperature in 
which our model might be qualitatively compared with QCD. 
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